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Abstract. The Hubble constant Hq describes not only the expansion of local space at redshift 
z ~ 0, but is also a fundamental parameter determining the evolution of the universe. Recent 
measurements of Hq anchored on Cepheid observations have reached a precision of several 
percent. However, this problem is so important that confirmation from several methods is needed 
to better constrain Ho and, with it, dark energy and the curvature of space. A particularly direct 
method involves the determination of distances to local galaxies far enough to be part of the 
Hubble flow through water vapor (H2O) masers orbiting nuclear supermassive black holes. The 
goal of this article is to describe the relevance of Ho with respect to fundamental cosmological 
questions and to summarize recent progress of the "Megamaser Cosmology Project" (MCP) 
related to the Hubble constant. 
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1. Cosmological Background 

For 85 years, it has been known that our universe is expanding (Lemaitre 1927). This 
expansion was beheved to slow down in time because of gravitational attraction. However, 
based on observations of luminous standard candles (Type la supernoave) Riess et al. 
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(1998) and Perlmutter et al. (1999) suggested instead accelerated expansion, turning 
cosmology upside down and winning the most recent Nobel Prize in physics. More than 
a decade after this discovery, accelerated expansion is well established. A de-acceleration 
during the initial few billion years after the "Big Bang" , when densities of matter and 
radiation were much higher than today, is followed by accelerated expansion. The cause of 
the accelerated expansion is so far unknown and is described by the term "Dark Energy" . 
Following the standard model, it should account for the majority of the energy density 
of the universe and retards the formation of large scale structure. Understanding dark 
energy may be the most important problem existing in physics today. 

Dark Energy dominates the energy budget, accelerates the expansion of the universe, 
and affects large scale structure. What is its nature? There are three classes of potential 
explanations: (1) a cosmological constant, which has been proposed already in the early 
days of general relativity (Einstein 1917) as a kind of repulsive gravity, (2) a scalar field, 
somewhat analogous to that proposed to explain inflation at a much earlier time (e.g., 
Wetterich 1988; Ratra & Peebles 1988), and (3) modified gravity (e.g., Tsujikawa 2010), 
which will not be considered here. 

Assuming that the universe is homogeneous and isotropic (as approximately suggested 
by the large scale matter distribution and the 3K microwave background), the space-time 
metric can be written in the following form 

ds"^ = dt'^ -a'^it) X [dr'^/{l-kr^)+rde^ +r^sm'^e d(l)^], (1.1) 

with t and a oc (1 + z)~^ being time and cosmic scale factor, r, 6, and denoting 
comoving spatial coordinates, and k representing the curvature of 3-dimensional space. 
The field equations of general relativity, applied to this Priedmann-Robertson- Walker 
metric, lead to the so-called Friedmann equations, 

„ fdV SttG , c^ A 

and 

a 47rG , „ , A „, 
a =-^(^ + ^^^+3- ^'-'^ 
H is the Hubble parameter for a given time {Hq stands for redshift z = 0), p is the 

density of matter and radiation, p denotes pressure, and A represents the traditional 
cosmological constant, which (like dark energy) can be subsumed into the density and 
pressure term, 

a 47rG v^, „ , 

- = -^ E(a + 3pO. (1.4) 

i 

Defining w = p/p, the different (in part putative) components yield: 



Matter: w = 
Radiation: w = 1/3 
"Quintessential" scalar field: -1 < w < -1/3 
Cosmological constant: w = -1 
Phantom energy: w < -1 



The resulting values of w determine the normalized acceleration d/a. For gravity (mat- 
ter) we obtain the expected negative value. This also holds for radiation, which has dom- 
inated, according to the standard model, at redshifts JjlO^. The {p + 3p) term in Eq. 1.4 
directly infers that w values smaller than -1/3 are required for accelerated expansion. 




Figure 1. Constraints from the cosmic microwave background (Spergel et al. 2007). The different 
colors correspond to different values of the Hubble constant Ho. is the energy density of 
baryonic and dark matter in units of the critical density, Qa is the corresponding parameter for 
dark energy. For a "fiat" universe, fim + S7a ~ 1. The data are consistent with a wide range of 
Ho values. 



which is therefore the possible w range for dark energy, w and H, and thus also Hq, are 
obviously related, emphasizing the cosmological importance of the Hubble constant. 

Fig. 1 shows the range of not necessarily flat cosmological cold dark matter (CDM) 
models consistent with Wilkinson Microwave Anisotropy Probe (WMAP) data (Spergel 
et al. 2007). Assuming that the universe is flat would provide a rather accurate solution, 
but is it really flat as suggested by inflationary models? The present value of the curvature 
radius, Rq, is related to Hq and = + Ha (see Fig. 1 for definitions) by 

Ro - (a/2) X |fcri/2 ^ (c/2)x H^^ x (no-l)-^^^. (f.5) 

c is the speed of light and (c/2) x Hq^ ^ 2.f Gpc is known as the Hubble radius. 
Obviously, with our present precision to determine flo, curvature radii as small 10 Gpc 
cannot be excluded and we are still far from being able to state that the universe is 
truly flat. All this implies that the cosmic microwave background alone does not provide 
stringent limits. This is not unexpected, since the CMB provides information at a single 
early epoch, when dark energy did not play a role. Its importance can only be deduced 
by a combination of CMB observations with data from the much younger universe. 

There are several lines of observational activity to constrain the wide range of models 
permitted by the CMB: These include (1) Type la supernovae, the standard candles, 
where luminosity distances and redshifts can be compared; (2) galaxy clusters, where 
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redshift independent ratios between baryonic and dark matter masses can only be ob- 
tained with a small subset of possible cosmologies; (3) gravitational lensing or cosmic 
shear, where the dark matter distribution can be determined, isolating dark energy; and 
(4) baryon acoustic oscillations (BAO), which left an imprint on the cosmic microwave 
background, which is clearly seen in CMB power spectra (e.g., Spergel et al. 2003, 2007; 
Komatsu et al. 2011). This imprint is still visible affecting structure at moderate red- 
shifts and providing a local size scale of order 140 Mpc. All these tracers are observed 
at significant redshifts. However, it is at redshift zero where dark energy is most domi- 
nant, and only here its energy density is significantly higher than that of baryonic and 
dark matter combined. Thus it is the Hubble constant, providing a measure of the local 
universe, which provides the longest lever arm with respect to the CMB to measure the 
effects of dark energy (Hu 2005). 

2. Constraining Hq 

Assuming a A cold dark matter universe and excluding curvature, exotic neutrino or 
specific early universe physics, Komatsu et al. (2011) derive from WMAP data Hq = 
71.0 ± 2.5kms^^ Mpc~^. While studies based on gravitational lens time delays (e.g., 
Treu & Koopmans 2002; Cardone et al. 2002) and on X-ray and Sunyacv-Zel'dovich data 
of galaxy clusters (e.g., Bonamente et al. 2006) have been used to constrain Hq, such 
deductions from rcdshiftcd objects depend on the chosen cosmological model and are no 
substitute for a measurement of Hq in the local universe. 

With the HST key project to measure the Hubble constant, Frecdman ct al. (2001) 
obtained from Cephcids in nearby galaxies Hq = 72 ± 3^ ± 7s kms^^ Mpc^^, estimat- 
ing both random and systematic errors. This was based on the extragalactic distance 
ladder using Cepheid variable stars in the Large Magellanic Cloud to calibrate the mea- 
surements. However, the LMC has a low mctallicity and significant depth. A potential 
dependence of Cepheid luminosities on mctallicity and uncertainties in the distances to 
individual stars complicate the calibration. This was highlighted by Sandage et al. (2006), 
who used similar methods but a different mctallicity correction to obtain Hq = 62 ± 1.3r 
± 5.0s kms~^ Mpc~^ (Fig. 2). While the Hq value by Freedman et al. (2001) is consistent 
with a flat universe, the Sandage et al. (2006) result challenges it. 

To date the most ambitious program to determine Hq is that led by A.G. Riess (see also 
Freedman & Madore 2010 for a recent review). Based on three anchors, (1) the parallax 
determinations of Galactic Cepheids, (2) Cepheides in the LMC witli distances deduced 
from eclipsing binaries, (3) the distance to NGC4258 (see Sect. 3) and with new HST 
(Hubble Space Telescope) data from galaxies with Cepheids and Type la supernovae, 
Riess et al. (2011) derive Hq — 73.8 ± 2.4kms~'^ Mpc"'^. This estimate also makes use 
of Cepheids measured in the near infrared, which helps to reduce both systematic and 
random errors with respect to optical observations. Note that the distance to NGC 4258 
has been slightly revised in the meantime (E.M.L. Humphreys, priv. comm.). 

While all these measurements are highly encouraging and pave the way to a more pre- 
cise knowledge of our universe, a totally independent measure of Hq is essential to either 
confirm the above cited results or to hint at problems that may have been overlooked so 
far. 

3. The Megamaser Cosmology Project (MCP) 

The MCP is an NRAO (National Radio Astronomy Observatory) key project to de- 
termine Hq by measuring geometric distances with an accuracy of ~10% to ^lO galaxies 
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Figure 2. WMAP (Wilkinson Microwave Anisotropy Probe) la and 2a likelihood surfaces for 
f2m (see Fig. 1) and w (Eq. 1.4), with priors on Ho (upper panels: 72kms~^ Mpc~^; lower panels: 
62kms~^ Mpc). Right versus left panels (solid lines) demonstrate the improvements gained by 
reducing the uncertainty in Ho from 10% to 3%. Dotted lines: "wcdm + no perturbations" model 
from Spergel et al. (2007); solid lines: the same, but with constraints from Hq incorporated. 



in the local Hubble flow. Following the experience gained by studying the prototypical 
source, NGC 4258 (Miyoshi et al. 1995; Herrnstein et al. 1999), it includes (1) a GBT 
(Green Bank Telescope) survey to identify suitable circumnuclear 22 GHz H2O maser 
disks (see Fig. 3 for a spectrum), (2) direct imaging of these sub-pc disks, using the 
VLB A (Very Long Baseline Array), the GBT, and, for northern sources, also the Eff els- 
berg telescope, (3) GBT monitoring to measure accelerations of the spectral components, 
and (4) model calculations to simulate; the maser disk dynamics. So far published arti- 
cles include Reid et al. (2009), Braatz et al. (2010), Greene et al. (2010), and Kuo et al. 
(2011). 

The 22 GHz H2O "megamasers" , luminous masers associated with active galactic nu- 
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Figure 3. One of the newly found maser disks with systemic (green-yeUow) and high velocity 
(violet-blue and orange-red) components. 



clei, are mostly found in Seyfert 2 and LINER (Low Ionization Nuclear Emission Line 
Region) galaxies with high column densities (Braatz et al. 1997; Zhang et al. 2006; Made- 
jski et al. 2006; Greenhill et al. 2008), relatively high optical luminosity, velocity disper- 
sion, and [Oiii]A5007 luminosity (Zhu et al. 2011) as well as relatively strong Fe Ka lines 
in those sources which are Compton thick (Zhang et al. 2010). Low X-ray/[Oiv]A25890 
ratios (RamoUa et al. 2011), and high nuclear radio continuum luminosities (Zhang et 
al. 2012) with respect to II2O undetected galaxies are also statistically obtained. 

The H2O masers reside in thin, edge-on gaseous annuli. Emission near the systemic 
velocity of the parent galaxy originates from the near side of the disk and red- and blue- 
shifted satellite lines come from the two tangent points (see Figs. 3-5). Assuming an 
ideal circular, warpless thin disk, seen perfectly edge-on, and in Keplerian motion, the 
mass Magn enclosed by the disk is 



Magn = 1.12 X 



kms 



R 



mas 



D 
Mpc 



M 



0j 



(3.1) 



with T^-ot denoting the rotation velocity at angular radius R and D representing the 
distance. Very long baseline interferometry maps allow us to directly measure VJ-ot for 
various values of R. From the Keplerian rotation curve we then obtain the constant 



Ci 



Vrot 

kms" 



R 



1/2 



(3.2) 



The velocity gradient of the systemic features as a function of impact parameter provides 
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Figure 4. Left panel: H2O image of tlie maser disk in UGC 3789 (from Reid et al. 2009). 
Tlie insert presents a magnification of the systemic features. Right panel: Radial velocity versus 
impact parameter (also Reid et al. 2009). For the blue- and red-shifted high velocity components, 
Keplerian r~^^^ rotation curves are also displayed. 



another constant, 



C2 = 



kms 



R 



(3.3) 



C1/C2 = then gives the angular radius Rg of the systemic features as viewed from 
a direction in the plane of the disk, but perpendicular to the line of sight. The total 
distance to the galaxy is then determined by the centripetal acceleration 



dVs/dt 



rot 



(3.4) 



with the index "s" denoting the systemic maser components. With dVs/dt being mea- 
sured, the linear scale rg can be compared with the angular scale Rs to provide the 
preliminary distance estimate. 

While such a procedure leads to a rough first estimate, disks may be warped (see 
Miyoshi et al. 1995 for the first such case, NGC 4258), orbits may be eccentric, and 
inclinations may not equal 90°. To model the circumnuclear disks as seen in H2O as 
detailed as possible, a Bayesian fitting procedure has been developed (M.J. Reid), using 
a Markov Chain Monte Carlo approach. A Metropolis Hastings algorithm is applied to 
choose successive trial parameters covering the parameter space. Fig. 6 displays such a 
simulation for NGC 6264. These very preliminary simulations indicate low eccentricities 
(e < 0.1). 

For UGC 3789, Fig. 5 indicates two groups of systemic maser components, one with 
a higher acceleration than the other. Unlike in NGC 4258, the systemic features do not 
arise from a single ring segment with specific galactocentric radius. The distance to the 
galaxy can be derived separately for the two rings. These yield Di = 50.2 ± 7.7 Mpc 
and D2 = 48.1 ± 17.4 Mpc. The weighted mean is -Dugc3789 = 49.9 ± 7.0 Mpc (14%). 
With a peculiar radial velocity relative to the cosmic microwave background of -151 
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Figure 5. Acceleration of individual maser features in UGC 3789 
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Figure 6. Bayesian fitting results of the maser disk in NGC 6264. Left panel: Face-on View 
onto the disk with systemic (green), approaching (blue) and receding (red) components. Right 
panel: Model of the warped disk viewed approximately edge-on as observed from Earth. 

± 163kms~^ and Vcmb = Vlsr + 60knis~^, the relativistic, recessional flow velocity 
becomes 3481 ± 163kms~^. This yields with standard ACDM parameters Hq = 69 ± 
llMpc-i (Braatz et al. 2010). 

For NGC 6264, a source at a distance of 150 Mpc, a similar analysis results in = 
65.8 ± 7.2kms~^ (Kuo 2011). Combining both sources, our present best estimate for the 
Hubble constant becomes Hq = 67 ± 6kms~"'^. This preliminary result is so far consistent 
with all previously (Sect. 2) mentioned values. 



4. Prospects 

So far, two sources yielded publishable results. Increasing this number to 10 and ac- 
counting for the fact that the targets are located in different parts of the sky, the Icr 
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Figure 7. The maser disk of NGC 4258 (Argon et al. 2007), extended east-west, and below tlie 
maser dislc of NGC6323, extended approximately north-south, on the same angular scale. The 
given beam shows a typical synthesized beam for a high-declination target at 22 GHz using the 
VLBA, the GBT, and Effelsberg. The sketch emphasizes the great progress achieved in recent 
years when mapping H2O maser emission in distant sources. 

error of 6kms~^Mpc~^ obtained so far should decrease by a factor of (2/10)"'^/^ to 
~2.7kms~^ Mpc~^ or 4%. Longer monitoring and more interferometric maps can reduce 
this uncertainty further. To demonstrate the degree of sensitivity required for these mea- 
surements, Fig. 7 shows the size of the prototypical nuclear disk in NGC 4258 on the same 
angular scale as the disk toward NGC 6323. While NGC 4258 is with V ~ 500 km s'^ not 
yet in the Hubble flow and therefore not useful for a direct estimate (its maser lines 
are nevertheless essential to calibrate the distance scale defined by Cepheids), NGC 6323 
has the potential to probe the distance scale with its recesseional velocity of almost 
7800kms-i. Toward NGC 6323, NGC 1194, NGC 2273, and Mrk 1419 the maser disks 
have also been mapped, demonstrating that the technique to derive distances, first tried 
out on NGC 4258, can also be used for much more distant galaxies. 

Figs. 3, 5, and 7 directly demonstrate the importance of sensitivity. While in UGC 3789 
many features have flux densities of ~5 mJy or higher, which can be readily analyzed, 
more distant sources reveal a plethora of components below this critical level. Getting 
these components as well would greatly facilitate any analysis. Thus the inclusion of a 
phased Jansky VLA (Very Large Array) is highly desirable. The completion of the Sar- 
dinia telescope may also help in the foreseeable future. Furthermore, as Fig. 7 indicates, 
angular resolution is another essential point. While an SKA- high would guarantee ex- 
treme sensitivity, the small angular extent of the H2O maser disks requires a world-wide 
array, with space- VLBI providing another significant improvement. 

Aside of the maser sources mentioned above, new targets have been detected, which 
look promising when analyzing their single dish spectra. While it remains to be seen how 
useful they will be for detailed mapping, it is worth mentioning that so far no systemic 
feature has been detected that shows a secular drift to the blue side. Either the nuclei 
are opaque at 22 GHz or the radiation is so highly beamed that maser photons from the 
backside of the disks have no chance to reach us. 
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